Dynamics of the Destruction and Rebuilding of a Dipole Gap in Glasses 
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After a strong electric bias field was applied to a glass sample at temperatures in the millikelvin 
range its AC-dielectric constant increases and then decays logarithmically with time. For the 
polyester glass mylar we have observed the relaxation of the dielectric constant back to its initial 
value for several temperatures and histories of the bias field. Starting from the dipole gap theory 
we have developed a model suggesting that the change of the dielectric constant after transient 
application of a bias field is only partly due to relaxational processes. In addition, non-adiabatic 
driving of tunneling states (TSs) by applied electric fields causes long lasting changes in the dielec- 
tric constant. Moreover, our observations indicate that at temperatures below 50 mK the relaxation 
of TSs is caused primarily by interactions between TSs. 

PACS numbers: 61.43.Fs, 77.22.-d, 05.70.Ln 



At low temperatures the internal energy of glasses is 
vastly elevated by atomic tunneling states (TSs) and 
most properties of amorphous solids are strongly influ- 
enced by these additional degrees of freedom [[j] [|. The 
phenomenological 'tunneling model' describes the ther- 
mal, elastic and dielectric properties of glasses at low 
temperatures successfully || Q. Starting from the cen- 
tral assumption that the potential minima of groups of 
atoms are not well defined, it describes such a configu- 
ration by a 'particle' moving in a double well potential. 
The two identical and harmonic wells may differ in depth 
by the asymmetry energy A. At low temperatures only 
the energy ground states of the two wells are occupied 
and the TSs oscillate by quantum mechanical tunneling 
between them. The energy splitting of these two-level 
systems is given by Eq — y/ A% + A 2 , where Ao is the 
tunneling splitting. The parameters Ao and A are widely 
distributed due to the randomness of the glassy structure. 

Soon after the introduction of the tunneling model 
in 1972 it was shown that the mutual interaction of 
TSs causes spectral diffusion, while interaction that in- 
volves energy exchange between TSs was long believed 
to be unimportant. However, more recent experiments 
revealed deviations from the tunneling model and indi- 
cate that for many properties the interaction between 
TSs gains importance with decreasing temperature || |) . 

Here we will concentrate on the effects of electric bias 
fields on the low temperature AC-dielectric constant of 
amorphous materials. A sudden large change of the elec- 
tric bias field causes the relative dielectric constant Se/e 
of a glass sample, that was originally in thermal equilib- 
rium, to immediately change to a larger value and then 
decay back to its initial value logarithmically for several 
decades in time. In addition, a relatively slow sweep 
of the bias field reveals a minimum of Se/e at the ini- 
tial bias field, in our case at zero field ||. A. Burin ex- 
plained these findings within a 'dipole gap' model, which 
considers the weak dipole-dipole interaction Jy = Uo/r§ 
between TSs M . Although the interaction is weak, some 



TSs fulfill the condition Jy > (^oi+^Oj) f° r strongly cou- 
pled pairs. Within the dipole gap model, such strongly 
coupled TSs form clusters that do not contribute fully 
to the dielectric constant because their individual dipole 
moments are pinned relative to each other. A change of 
the bias field alters the distribution of asymmetry ener- 
gies and leads to the formation of new clusters of strongly 
coupled TSs. As long as the new clusters remain out of 
thermal equilibrium, their TSs contribute fully to the 
dielectric constant Q. This explains the observed mini- 
mum of Se/e at the initial bias field. In addition mutual 
interaction of TSs leads to a new relaxation channel || 
due to resonant energy exchange between pairs of TSs 
resulting in the overall relaxation rate 
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consisting of the well known one-phonon relaxation rate 
(first term r.h.s.) and the contribution due to inter- 
action. The material constants 70 || and ao S de- 
pend on the coupling of the TSs to phonons and be- 
tween each other respectively. Note that the dynamics 
of TSs that generate any change in the dielectric con- 
stant is dominated by thermally active TSs that fulfill 
E = y/Al + (A + p • F) 2 ~ k B T producing an addi- 
tional temperature dependence in (|l]) implicit in the en- 
ergy splittings E. Moreover, after changing the bias field 
F by more then p • 5F ~ k B T, where p is the dipole mo- 
ment of the TSs, a different set of TSs fulfills E ~ k B T. 
Therefore in past experiments in which a large bias field 
was suddenly applied to study the dynamics of strongly 
coupled clusters the formation of a new dipole gap con- 
sisting of different clusters than those formed at zero field 
was observed. The interpretations of these measurements 
relied on assumptions such as a flat energy distribution 
of TSs for E < p • F = $. 

All our measurements were made at zero bias field. 
Thus we always measured the same set of TSs. However, 
we applied a large electric bias field $ 3> k B T for the 
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waiting time t w to destroy the zero field dipole gap and 
then observed the restoration of the same dipole gap by 
measuring Se/e in zero field. In addition, we took care 
to linearly sweep the bias field up and down in the same 
sweep time t s . 

As long as the bias field is applied the TSs are relaxing 
towards a new equilibrium, thereby breaking clusters of 
TSs that were initially strongly coupled. Back at zero 
bias field these clusters rearrange. The TS's relaxation 
times depend according to (|l|) on their asymmetry ener- 
gies and are broadly distributed. For $ » k B T initially 
thermal TSs with E ~ k B T have E ~ $ at the bias field 
F. We introduce the decay time tq at which we expect 
the dielectric constant to again reach its initial value, i.e. 
e(t) — e(oo) before the bias field was applied 

to = t{E = k B T) _ fk B T\ 2 ( 7o $ 3 + a k B T) 
t w t{E = ~ \ J 2.2 7o (fc B r)3 + a k B T ■ 

(2) 

For this we assumed that both the destruction and 
restoration of the dipole gap are solely caused by relax- 
ational processes, thus that tq is given by the ratio of the 
relaxation times of TSs without and with field times the 
waiting time. Following the dipole gap theory we obtain 
for the change of the dielectric constant, after a bias field 
F was applied for a time t w , 

Since in ^ rclaxational processes are approximated by 
Heaviside step functions we expect deviations from loga- 
rithmic decays at the start and the end of the decay. 

Our experiments were performed on 15/xm thick my- 
lar film (amorphous polyester) at a frequency of 10 kHz 
at temperatures well below the minimum of Se/e at 
T m in(10kHz) ~ 50 mK, ensuring that Se/e consists en- 
tirely of its resonant part. Gold capacitor electrodes were 
evaporated onto the sample. The dielectric constant was 
measured with a hand made bridge by comparison to a 
vacuum reference capacitor, that was placed on the one- 
kelvin pot of the dilution refrigerator. The AC-field was 
always kept within its linear response region. 

Fig. [j] shows the evolution of Se/e after a bias field 
of F = 6.7 MV/m was applied for two different waiting 
times of several hours (t w ^> t s = 1 s) at each of the 
temperatures T — 10 mK and T — 20 mK. The zero 
point of time (t = 0) is chosen to be at the completion 
of the bias field sweep and Se/e = is defined as its 
initial value at zero field. The amplitude of Se/e ~ 10~ 5 
increases for longer waiting times and decreasing temper- 
ature as expected from (||). The dielectric constant de- 
cays approximately logarithmically over several decades 
in time justifying the simplification of describing relax- 
ational processes by Heavyside Step Functions. Taking 
only the one-phonon process into account («o = 0) we 
would expect decay times of tq ~ 73i w for T = 10 mK 
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FIG. 1: Relative dielectric constant of Mylar after a bias field 
F was applied at two temperatures for very long waiting times 
iw S> t s . All lines are theory curves after (see text). The 
shematic indicates the application of the bias field and the 
choice of t = 0. 



and To ~ 37t w for T = 20 mK according to (Q). In con- 
trast, our measurements suggest decay times in the order 
tq ~ i w indicating the importance of interaction medi- 
ated relaxation ]lQ| . Furthermore our experiments indi- 
cate that the restoration of the dipole gap takes longer 
for the higher temperature. According to ([!]) it is the 
interaction mediated relaxation process that results in 
an increasing relaxation time for the contributing TSs 
(Eo ~ k B T) with increasing temperature. Therefore our 
data unambiguously show that interaction mediated re- 
laxation dominates in this temperature range. 

One of the solid lines in Fig. [l] is a linear fit for 
T = 10 mK and t w = 18000 s (filled circles). By com- 
parison with (^|) we find the fit-parameters r /t w ~ 1.02 
and PoIIq ~ 4.1 • 10 -4 , the latter lying whithin typical 
estimations for other glasses Q. The prefactor in (||) was 
previously obtained from the temperature variation of 
Se/e well below T min as 2P ^ ~ 1.35 • 10~ 4 §, [ll| and 
the dipole moment of the TSs in mylar was estimated to 
be p ~ 1.2 Debye as discussed later (see fig. ||). Using (|l|) 
we estimate the temperature at which the interaction me- 
diated relaxation time equals the one-phonon relaxation 
time as 55 mK < T c < 95 mK for 3 > T /t w > 1. Be- 
low T c the interaction mediated relaxation dominates. 
The other three lines in fig. [j] are calculated from (||) us- 
ing the same set of parameters. Our simplifying model 
agrees remarkably well with the data besides deviations 
from logarithmic decays particularly at the higher tem- 
perature. 

In fig. |2| the restoration of the dipole gap at the temper- 
ature of T = 10 mK is displayed for 10 different waiting 
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FIG. 2: Relative dielectric constant of Mylar after a bias field 
F was applied for waiting times i w as marked left of each 
curve, while t s = 1 s for all data. The solid lines are theory 
curves after (|) using the same parameters as in fig. |l| 



times between 0.1 s and 5 hours and otherwise the same 
conditions as for the data in fig. |l|. The straight lines in 
fig. H are calculated from (||) for the three longest wait- 
ing times using the same parameters as for the theory 
curves in fig. |} Our model which considers relaxation 
predicts parallel decay curves shifted along the logarith- 
mic time axis proportional to ln(i w ). In experiments we 
observe this behaviour only for our longest waiting times 
(fig. 0). As t w gets shorter Se/e increasingly decays slower 
than expected, until for t w < t s (our two shortest waiting 
times in fig. |^) the waiting time dependence disappears 
completely. To further investigate this puzzling effect 
we show in Fig. ^ the decay curves for 3 different sweep 
times at t w = 30 s (solid points). As with longer t w , Se/e 
increases for longer t s . For comparison we included the 
data points from Fig. | for t s = 1 s and t w — 100 s (open 
circles). Surprisingly their course is almost identical to 
the curve for t s = 10 s and t w = 30 s with a much smaller- 
waiting time even if one were to make the obvious over- 
estimation t w + 2t s — > t w to include relaxation while the 
bias field is being changed. 

At this point we can safely conclude that at waiting 
times shorter than t w ~ 3000 s an additional process 
gains importance that is not relaxational (due to waiting 
time independence) and instead depends on the sweep 
rate of the bias field. We propose that the changing 
electric (AC and DC) field drives TSs non-adiabatically 
during the bias field sweep, and that this process is re- 
sponsible for the new behavior. 

Let us consider the time dependent quantum mechan- 
ics of the TSs while the bias field is being changed. If 
the field sweep is 'slow' the TSs are driven adiabatically 
and their wave functions always adjust to the momentary 
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FIG. 3: Relative dielectric constant of Mylar after a bias field 
F was applied for different t B and t w . 



field. In the opposite case of a 'fast' sweep the TSs are 
driven non-adiabatically and their wave functions stay 
unchanged. These TSs remain in their original quantum 
states. Most TSs are asymmetric and therefore localized. 
In the following we consider non-adiabatically driven TSs 
having asymmetry energies that decrease to zero and suc- 
cessively increase again with changing field. For these 
TSs the original ground (excited) state at zero field con- 
verts into the excited (ground) state at maximum field. 
As the result the occupation numbers of these TSs are in- 
verted. However, as the field is changed back to zero the 
non-adiabatically driven TSs are in their original state 
having been inverted twice. The Landau-Zenner crite- 
rion |12f| gives the probability 



fli(Ao) = e 



with A 0c = TTy/^h/U (4) 



with Plank's constant h that a system is driven non- 
adiabatically near A + $ = 0. Only TSs with Ao ~ Ao c 
have a finite probability to have been inverted exactly 
once 1 13 . The non-adiabatically driven TSs have very 
long relaxation times r(A 0c ,-E'o) which are according to 
(0) and (jj) proportional tor~ t s /&. Neglecting the re- 
laxational process leading to (||) we expect the decay time 
at which the dielectric constant reaches its initial value to 
be a function of t s /& as well. However, the predicted de- 
cay would be non-logarithmic resembling a smoothened 
step occuring at the time t ~ r(Ao c , Eq — k^T). 

An AC measuring field forces the dipole moment of a 
TS to oscillate with the measuring frequency lj as long 
as the quantum mechanical time evolution of the TS is 
able to follow the field. However, this description fails 
when the energy splitting E of a TS is smaller than a 
critical value E\c = P ' Fac below which the AC- 
measuring field F^c drives the system non-adiabatically. 
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TSs with Ao < Eac are driven nonadiabatically by 
the AC field for the time St* = t s Fac I -Pbc m which 
E < Eac ■ An initially localised TS will oscillate between 
the two wells during St* with a renormalizcd tunneling 
frequency A /h @. Accordingly, TSs with hA~ 1 > St* 
stay localised in the original state whereas TSs with 
/iAq 1 <C St* oscillate many times so that only half of 
them end up in their initial state. The AC-driving ef- 
fect predicts a logarithmic decay with the decay time 
T ((Si*~i B /$) p. 
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A combination of DC- and AC-driving explains the 
sweep time dependence visible in fig. || and || qualita- 
tively. In our experiments -Fac — 9.7 kV/m, resulting in 
Eac ~ A 0c for F DC ^ 6.7 MV/m and t s = 1 s. At longer 
t s we expect AC-driving to increasingly become impor- 
tant at the expense of the DC-driving effect. The data 
in fig. H indicate a crossover from non-logarithmic to log- 
arithmic behavior at t s ~ 10 s. Consistently we did not 
observe a dependence on Fac for t s < 10 s (unpublished 
data) but have not searched for that effect for t s > 10 s. 
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FIG. 4: Relative dielectric constant of Mylar after a bias field 
F was applied for the same sweep rate, the same waiting time 
but different bias field strengths and sweep times. In the inset 
the data are normalized with respect to their bias field for the 
dipole moment of p — 1.2 D. 

Our model including non-adiabatic driving predicts 
t = To(is/$Dc)- in n g- H we pl°t t w0 data sets with 
the same 3>/i s = 19 fcemK/s but different sweep times 
and bias fields. The decay times of both data curves are 
the same as expected. We stress that relaxational pro- 
cesses during £ w cause a decay time To ~ t w (fig. |l|), but 
in fig. [|to ^> t w and thus the decay is due to the driving 
effect. In our model this can be understood by recog- 
nizing that the finite waiting time necessarily implies a 
smallest Ao for TSs that can possibly decay while the 
bias field is applied. This smallest Ao is larger than Ao c 
unless t w ^> t s . Regardless of which effect breaks the 



strongly coupled clusters of TSs, we always find the tem- 
perature dependence of the decay of the dielectric con- 
stant propotional to [ln($/fcBT)] 2 . This is a consequence 
of the central assumption in the dipole gap model that 
only TSs with energy splitting T < Eq < $ contribute to 
changes in the dielectric constant ||]. Using this depen- 
dence we normalized the two data sets as shown in the 
inset in fig. ^. Very good agreement was reached with 
the dipole moment p ~ 1.2 Debye (that was used for the 
theory curves in fig. |l| and ||). 

We have investigated the dynamics of the destruction 
and the restoration of the zero field dipole gap by appli- 
cation of strong electric bias fields. Only our data curves 
at very long waiting times t w 3> t s can be explained by 
solely relaxational processes. Our investigation suggests 
non-adiabatic driving of the TSs as a possible reason for 
the destruction of the dipole gap for shorter t w . Further- 
more we have shown that at temperatures below 50 mK 
the relaxation of TSs is dominately by a relaxation mech- 
anism due to interactions between TSs 
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